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We present a new QCD event generator for hadron collider which can calculate one-, two- and 
three-jet cross sections at next-to-leading order accuracy. In this letter we study the transverse 
energy spectrum of three-jet hadronic events using the k± algorithm. We show that the next-to- 
leading order correction significantly reduces the renormalization and factorization scale dependence 
of the three-jet cross section. 
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The latest version of the experiment at Tevatron and 
the future collider LHC will provide precise data so that 
not only inclusive measurements can be used to study the 
physics of hadronic final state. The studies of the event 
shapes and multi-jet event can be important projects. 

One of the important theoretical tools in the analysis 
of hadronic final states is perturbative Quantum Chro- 
modynamics (QCD). In order to make quantitative pre- 
dictions in perturbative QCD, it is essential to perform 
the computations (at least) at the next-to-leading order 
(NLO) accuracy. In hadron collision the most easily cal- 
culated one- and two-jet cross sections have so far been 
calculated at NLO level |0, ^. At next-to-leading level 
some three-jet observables were calculated by Giele and 
Kilgore ^, ^ and by Trocsanyi Q. In this letter we 
present a new NLO event generator for calculating jet 
observables in hadron-hadron collision. We compute the 
three-jet cross sections using the k±_ algorithm Q to 
resolve jets in the final state. With our Monte Carlo 
program one can compute the NLO cross section for any 
other infrared safe one-, two- and three-jet observables. 
The presented distributions are given simply as illustra- 
tion. 

In the case of one-jet inclusive cross section in high 
Pt region the forthcoming experimental data requires 
the knowledge of the higher order corrections. Recently, 
progress has been made in calculating the two loop 2^2 
matrix elements H, ^ . These matrix element are need 
to set up a Monte Carlo program for calculating the next- 
to-next-to-leading order (NNLO) one- and two-jet cross 
section. However, the three-jet NLO calculation, com- 
prising the one-loop 2^3 and tree-level 2^4 processes 
is a necessary first step of this project. A numerically sta- 
ble and fast three-jet NLO program can provide enough 
hope that it might be possible to develop a numerically 
stable Monte Carlo program for calculating one- and two- 
jet cross sections at NNLO level. 

In the last few years the theoretical developments 
make possible the next-to-leading order calculation for 
the three-jet quantities. There are several general meth- 
ods available for the cancellation of the infrared diver- 
gences that can be used for setting up a Monte Carlo 



program |jrT|, 12, |T^. In computing the NLO correction 
we use the dipole formalism of Catani-Seymour that 
we modified slightly in order to have a better control on 
the numerical calculation. The main idea is to cut the 
phase space of the dipole subtraction terms as introduced 
in Ref. 10, the details of how to this to apply for the 
case of hadron-hadron scattering will be given elsewhere. 

The advantages of using the dipole method are the fol- 
lowings: i) no approximation is made; ii) the exact phase 
space factorization allows full control over the efficient 
generation of the phase space; iii) neither the use of color 
ordered subamplitudes, nor symmetrization, nor partial 
fractioning of the matrix elements is required; iv) Lorentz 
invariance is maintained, therefore, the switch between 
various frames can be achieved by simply transforming 
the momenta; v) the use of crossing functions is avoided; 
vi) it can be implemented in an actual program in a fully 
process independent way. 

In this calculation we used the crossing symmetric tree- 
and one-loop level amplitudes. The parton subprocess 
-> 99999 II, ^ qqggg Q, ^ qqQQg and 
the subprocesses related to these have been computed to 
one-loop and ^ gggggg, qqgggg, qqQQgg 
p^ , pO| ^ and the crossed processes have been computed 
at tree level. 

We have checked numerically that in all soft and 
coUinear regions the difference of the real and subtrac- 
tion terms contain only integrable square-root singulari- 
ties. Furthermore, we have also checked that our results 
are independent of the parameter that controls the vol- 
ume of the cut dipole phase space, which ensures that 
indeed the same quantity has been subtracted from the 
real correction as added to the virtual one. 

Finally, to have a further check of the computation. 
Our method of implementing the dipole substraction 
terms allows for the construction of a process indepen- 
dent programming of QCD jet cross sections at the NLO 
accuracy. We use the same program structure, with triv- 
ial modifications, to compute one-, two- and three-jet 
cross sections. 

In order to check the base structure of the program we 
compare the our inclusive one-jet NLO result to predic- 
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FIG. 1: Comparison of the K factors of the one-jet inclusive 
cross section defined using the k± and for MRSD'_ parton 
densities obtained with Monte Carlo programs JETRAD and 
NLOJET+-I- (this work). The bands indicate the statistical 
error of the calculations. 

tion of the program JETRAD In this comparison we 
compare the one-jet inclusive cross section using the k± 
jet algorithm and MRSD'_ parton distribution function 
p2| . We find good agreement between the two program. 
The differences are within the statistical error as Fig. |l| 
shows. 

Historically, only the cone algorithm has been used to 
reconstruct the jet at hadron collider. In the three or 
higher jet calculation at NLO level the cone algorithm is 
not suitable because it has a lot of difficulties: an arbi- 
trary procedure must be implemented to split and merge 
the overlapping cones, and an ad- hoc parameter Rsep is 
required to accomodate the difference between the jet 
definitions at parton and detector level. To avoid this 
uncertainty we use the fc^ algorithm which has been de- 
veloped by several groups |^. Our implementation is 
based on the Ref . [|| . The algorithm starts with a list of 
the particles and the empty list of the jets. 

1. For each particle (pseudo-particle) i in the list, de- 
fine 

dz = Pt.i ■ (1) 
For each pair (i, j) of momenta (i ^ j), define 
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(2) 



^ is square of the 



where ARfj = (ry^ - rjj)^ + - 
angular separation which is expressed in the term 
of the pseudo-rapidity rji and the azimuth angle 0^ . 
£) is a free parameter. The usual choice of this 
parameter is D = I. 

2. Find the minimum of all the di and dij and label 
it dijiiii . 

3. If dmin is dij, remove particles (pseudo-particles) i 
and J from the list and replace them with a new. 
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FIG. 2: The perturbative prediction for the three-jet differ- 
ential cross section in the term of the transverse energy of the 
leading jet at Born level (light gray band) and next-to-leading 
order (dark gray band). The bands indicate the theoretical 
uncertainty due to the variation of the renormalization and 
factorization scales xr^f between 0.5 and 2. The solid line is 
the NLO result for the xa = xp = ^ choice of the scales. 



merged pseudo-particle given by the recombi- 
nation scheme. In this paper we use the E recom- 
bination scheme which is define the new pseudo- 
particle as the sum of the two particle 

P(ti) ^Pi+ P] ■ (3) 

4. If dmin is di, remove particle (pseudo-particle) i 
from the list of particle and add it to list of jets. 

5. If any particles remain, go to step 1. 

The algorithm produces a list of jets, each separated by 

Ai?y > D. 

Once the phase space integrations are carried out, we 
write the NLO jet cross section in the following form: 



njct 



^ J dl]adl]bfa/A{Va, fJ'%)fb/B{Vb, fJ-^) 
a.b 

[Pa,Pb,asi^lji), iir/Qh, fi^/Qn] , (4) 



where fi/H^Hit^'p) represents the patron distribution 
function of the incoming hadron defined at the factoriza- 
tion scale fip = xfQhi Va,b is the fraction of the proton 
momentum carried by the scattered partons Pa,b, Qh is 
the hard scale that characterizes the parton scattering 
which could be Et of the jet, jet mass of the event, etc 
and iiR = xrQh is the renormalization scale. 

Eq. (||) shows that using the dipole method one may 
either compute the full cross section at NLO accuracy 
including the convolution with the parton distribution 
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FIG. 3: The dependence of the three-jet cross section asjct 
on the renormahzation and factorization scales. 



functions, or simply the parton level cross section a at 
which can then be convoluted with the parton densities 
after the Monte Carlo integration. The latter procedure 
is the proper one if we are interested in measurement of 
the parton distribution functions (to avoid the recalcula- 
tion of the Monte Carlo integrals after each step of the 
fitting iteration). 

The three-jet cross sections presented here were calcu- 
lated for TEVATRON collider in proton-antiproton colli- 
sion at the center of mass energy = 1800 GeV. We re- 
strict the pseudo-rapidity range and the minimum trans- 
verse energy of the jets in laboratory frame to be 

- 4 < ?7jet < 4 , £;t > 50 GeV . (5) 

We choose the transverse energy of the leading jets 

(6) 



Qh — E^t'' 



as the hard scattering scale. 

In Fig. H, we plotted differential cross section in the 
term of the the transverse energy of leading jet convo- 
luted with the CTEQ5M1 parton distribution functions 
[ p3[ and using the two-loop formula for the strong cou- 
pling, 



as{M, 



^ as{Mz) Pi ln(w(^)) '^ 



/ N 1 as{Mz) 

win) = 1 - /3o — In 

Ztt \ 



27r Po w{p) 



(8) 



where as{Mz) = 0.118 and /3o = (HCa - iTRNf)/3, 
f3i = (17Ci - GCfTrNj - 10CATRNf)/3, with Nf = 
5 flavors. For the leading order results we used the 
CTEQ5L distributions and the one-loop ag {as{Mz) = 
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FIG. 4: The dependence of the three-jet differential cross 
section on the parameter D. The Rd means the ration of the 
differential cross sections for a given D and for D — 1. Upper 
figure shows the Born level result and the lower figure shows 
the NLO prediction. The error bars indicate the statistical 
error of the Monte Carlo calculation. 



0.127 and /?i in Eq. (0)). In this figure the theo- 
retical uncertainty of the three-jet cross section is shown. 
Over the wide range of the value the renormahzation and 
factorization scale (0.5 < xr^f < 2) this uncertainty is in 
the next-to-leading order result is much than in the Born 
level calculation. 

In Fig. H we study the scale dependences of the three- 
jet cross section. The strong dependence on the renor- 
mahzation scale observed at LO is significantly reduced. 
The factorization scale dependence is not significant at 
NLO and does not change much. Setting the two scales 
equal, fiR = = Ai, we can observe a wide plato peaking 
around xr — xf — 0.7. 

The inset figure in Fig. |^ shows the K factor (ratio of 
the three-jet cross section at NLO to that at LO accu- 
racy), indicating the relative size of the correction. We 
can see the size of the NLO correction is between 10 and 
25% for smaller values of the transverse energy and at the 
end of the spectra the size of the correction is almost zero. 
The error bars indicate the statistical error of the Monte 
Carlo calculation. Because of the strong logarithmic be- 
havior of the cross section the Monte Carlo calculation 
is very sensitive to the "missed binning" (when a huge 
positive comes from the real term and the corresponding 
huge negative weight form the subtraction term are filled 
in different histogram bins). 

In Fig. ^ we study the dependence of the differential 
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cros s section on parameter D. We plotted the ratios of 



the :ross section 



'dalf{E'4^-\) 



dE. 



(1) 



dE, 



(1) 



(9) 



for three different values of the parameter D {D = 
0.5, 1.5, 2). The parameter D controls the angular sep- 
aration in the jet algorithm procedure in Eq. ^. Chang- 
ing this parameter we expect more resolved jets (more 
three-jet events) with high transverse energy and less re- 
combination for smaller values of D and vice-versa. This 
behavior can be clearly observed from in Fig. ^. 

In this letter we presented a NLO computation of the 
three-jet cross section defined with the k± clustering al- 
gorithm in in hadron-hadron collision. Our results were 
obtained using a partonic Monte Carlo program that is 
suitable for implementing any detector cuts. We found 
that the NLO the correction is under 30% in the case 
of differential cross section but the K factor is sensitive 
to the allowed kinematic region. We demonstrated that 
the NLO corrections reduce the scale dependence signif- 
icantly. We also presented how the differential cross sec- 
tion depends on the angular separation parameter D used 
to define the jet. The same program can be used for com- 
puting the QCD radiative corrections to the (differential) 
cross section of any kind of one-, two-, or three-jet cross 
section or event-shape distribution in hadron-hadron col- 
lision. We compared the two-jet results obtained by our 
program to previous results and found agreement. 
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